Interplay between temperature and trap effects in one-dimensional lattice systems of 

bosonic particles 
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We investigate the interplay of temperature and trap effects in cold particle systems at their 
quantum critical regime, such as cold bosonic atoms in optical lattices at the transitions between 
Mott-insulator and superfluid phases. The theoretical framework is provided by the one-dimensional 
Bose-Hubbard model in the presence of an external trapping potential, and the trap-size scaling 
theory describing the large trap-size behavior at a quantum critical point. We present numerical 
results for the low-temperature behavior of the particle density and the density-density correlation 
function at the Mott transitions, and within the gapless superfluid phase. 

PACS numbers: 67.85.-d,05.30.Rt, 05.30.Jp 



Many-body phenomena in dilute atomic gases can be 
investigated in experiments of cold atoms in optical lat- 
tices created by laser-induced standing waves, such as 
Mott- Hubbard transitions in bosonic gases, see, e.g., 
Refs. [1-9]. An important feature of these experiments 
is the presence of a confining potential which traps the 
particles within a limited spatial region of the optical 
lattice. The theoretical framework [lfj is based on the 
Bose-Hubbard (BH) model [TO] 



(ij) 



22ni + 22v(ri)ni, 



(1) 



where (ij) is the set of nearest-neighbor sites, and m = 
b\bi is the particle density operator. We assume a power- 
law space dependence for the external potential which 
gives rise to the trap, i.e. 



V(r) = v p r p , 
with even p. The trap size is defined as 



(2) 



(3) 



Far from the origin the potential V(r) diverges, therefore 
(rii) vanishes and the particles are trapped. Experiments 
are usually set up with a harmonic potential, i.e. p = 2. 
In the following we set J = 1, so that I — l/v. 

The homogeneous BH model without trap undergoes 
continuous quantum transitions between superfluid and 
Mott-insulator phases, driven by the chemical potential 
fi. In one and two dimensions the length scale of the crit- 
ical modes diverges with the exponent v = 1/2 and the 
low-energy spectrum scales with the dynamic exponent 
z = 2 [ill II2I], The superfluid phase is gapless; in one- 
dimensional systems the low-energy modes are described 
by a two-dimensional conformal field theory. 

In the presence of a trapping potential, theoretical and 
experimental results have shown the coexistence of Mott- 
insulator and superfluid regions when varying the total 



occupancy of the lattice, see, e.g., Refs. @, E3, 
However, at fixed trap size, the system does not dev elop 
a critical behavior with a diverging length scale [HI, Il7j . 
Criticality can be recovered in the limit of large trap 
size, by simultaneously tuning the chemical potential to 
the critical values [i c of the homogenous system. This 
critical regime can be described in the framework of the 
trap-size scaling (TSS) theory [2^, [24|. At a quantum 
critical point, the large trap-size behavior of the free- 
energy density of a d-dimensional trapped particle system 
is expected to behave as [23|, |25[ 

F(ji, T, I, x) = r e(d+z) ^( j u/ e/l/ , Tl ez , xr e ), (4) 

where x is the distance from the middle of the trap, 
p, = fi — \x c and \i c is the critical value of the chemical 
potential. The trap exponent 9 at the Mott transitions 
is given by 



P 



P- 



(5) 



which determines how the length scale of the critical 
modes diverges with increasing trap size, i.e. £ ~ l e , 
at the critical point. Analogously, the particle density 
correlator, 



G(x,y) = (n x n y ) c = (n x n y ) - (n x )(n y ), 



(6) 



is expected to scale as 

G(x, y) = r 2de g(fll e ^, Tl ze , x/l , y/l e ). (7) 

Finite-size effects, when the trap is within a finite box 
of size L, can be taken into account by adding a further 
dependence on the scaling variable L/l B (26|. 

The TSS of the one-dimensional (ID) BH model has 
been investigated at zero temperature [25j |. TSS can 
be derived analytically in the low-density regime, at 
the superfluid-to-vacuum transition, within the spinless 
fcrmion representation of the hard-core limit. The trap- 
size dependence turns out to be more subtle in the other 



2 



critical regions, when the corresponding homogenous sys- 
tem has a nonzero filling /. In the presence of the trap- 
ping potential, the expectation value N of the particle- 
number operator N — n x is finite, and increases as 
N ~ I with increasing the trap size I keeping p fixed. 
Thus, since [N, Hbb] = 0, the lowest-energy states at 
fixed p, and in particular the ground state, change when 
iV jumps by 1 with increasing /, giving rise to an infinite 
number of level crossings [25[. At the Mott transitions 
with nonzero integer filling, this gives rise to a modu- 
lated TSS: the TSS of the observables is still controlled 
by the trap-size exponent 8, but it gets modulated by 
periodic functions of the trap size, requiring a minor re- 
vision of the TSS ansatz (j4]) and ([7]). For example, the 
particle-density correlation function at T — behaves as 
G(0,x) = l~ g(x/l e ,<p) where is a phase-like variable 
measuring the distance from the closest level crossing. 
Note that this phenomenon persists in the limit p — > oo , 
which corresponds to a homogeneous system with open 
boundary conditions. Modulations of the asymptotic 
power-law behavior is also found in the gapless super- 
fluid region, with additional multiscaling behaviors. 

In this paper we investigate the interplay of tempera- 
ture and trap effects in the ID BH model (JTJ at the transi- 
tions between Mott-insulator and superfluid phases. The 
results are analyzed and discussed in the theoretical TSS 
framework, which provides an effective description of the 
critical behaviors in the presence of the trap. In partic- 
ular, we address the issue whether the modulation phe- 
nomenon found at p < 1 and T — persists at finite 
temperature, or it is effectively averaged out due to the 
fact that the level crossings of the lowest states are ex- 
pected to give weaker effects at T > 0. This study is of 
experimental relevance, because quasi ID trapped par- 
ticle systems have been realized in optical lattices, see, 
e.g., Refs. [3,11 1. 

We present numerical results for the ID BH model, 
obtained by quantum Monte Carlo (QMC) simulations 
at fixed trap size I, based on the stochastic series expan- 
sion [H [Hf- The main features of the TSS at the Mott 
transitions are expected to be universal with respect to 
the on-site coupling U (provided U > 0), including the 
hard-core (HC) limit U oo. Therefore, we consider the 
ID HC BH model [29] , which is particularly convenient 
because it restricts the values of the site particle number 
to n x — 0,1, and it is expected to minimize the correc- 
tions to the universal TSS, as verified at T = [30|. In 
the absence of the trap, the ID HC BH model has three 
phases: the empty state for p > 1 with (n,) = 0, which 
may be seen as a particular n = Mott phase, a gapless 
superfluid phase for \p\ < 1, an d a n = 1 Mott phase 
for p < —1. See, e.g., Ref. 
two quantum transitions at p 



311 ] . Therefore, there are 
1 and p = —1, which 
share the same critical exponents v — 1/2 and z = 2. We 
present finite-temperature results for the particle density 
and the density-density correlator, at the transitions be- 
tween the Mott and superfluid phases, and within the 
gapless superfluid phase. Their temperature and trap- 



size dependences appear well described by the simplest 
TSS ansatz, such as that given in Eq. ([7|). The finite- 
temperature data do not show evidence of the modula- 
tion phenomenon, which appears averaged out, showing 
only a plain power-law TSS behavior. 

The paper is organized as follows. In Sec. [I] we study 
the temperature and trap-size dependences in the low- 
density regime, where we analytically derive the TSS de- 
scribing temperature and trap effects, and compare it 
with QMC simulations at p — 1. Sec. [H] presents QMC 
results at the n = 1 Mott transition for p = 2 and p — > oo. 
In Sec. |HT] we investigate the superfluid phase, present- 
ing QMC results at p = 0. In Sec. |IV] we study the 
TSS around the spatial point where the particle density 
rapidly vanishes when p < 1, which shows a peculiar 
scaling behavior effectively controlled by a linear exter- 
nal potential. Finally, we draw our conclusions in Sec.fVl 
In App. [A] we provide some details on the QMC simula- 
tions. 



I. THE LOW-DENSITY REGIME 

In order to determine the temperature and trap-size 
dependence in the low-density regime of the ID HC BH 
model, around p = 1, we consider the free fermion rep- 
resentation which can be derived by a Jordan- Wigner 
transformation, see, e.g., Ref. [31| . 



hij S-ij 



ij 



(8) 



with p = p — 1. In the fermion representation the 
Hamiltonian can be easily diagonalized by introducing 
new canonical fermionic variables r/k = Y]j 4>kiCii where 
4> satisfies the equation J2jhij<f>kj — ^k4>ki, obtaining 
H c = '^2i e ^ ! kVk r lk- The ground state contains all 77- 
fermions with LJk < 0. 

A nontrivial TSS limit [23j j around p = 1 is obtained 
by introducing the continuum function (f>k (x) = (j)k x , the 
rescaled quantities 



X = X/1\ p r =pl 2f> , n k =LO k l 



(9) 



where 9 is the trap exponent, cf. Eq. (|5|), and neglecting 
terms which are suppressed in the large-Z limit. This 
leads to the Schrodingcr-like equation 



1 d 2 



2 dX 2 



+ X p )ip k (X) = e k <p k (X), 



(10) 



where e k = £l k — p r and ip k (X) = 4> k {l e X). 

In the case of a harmonic potential, i.e. p — 2, we have 
6 = 1/2, and 

e k = n k ~ /i r = 2 1 / 2 (fc + l/2), fc>0, (11) 
MX) = ^f^y/l ex P (-XVV2), 
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FIG. 1: (Color online) The particle density (top) and the 
density-density correlator (bottom) in the presence of a har- 
monic potential, at fi — 1 and for some values of the trap size 
I, at r = Tl = 8 and r = 2. The full lines show the TSS 
functions given by Eqs. (|18|) and (|20[1 . 



where H^(x) are Hermite's polynomials. For p — > oo, 
Eq. (jTOJ) becomes equivalent to the Schrodinger equation 
of a free particle in a box of size L = 21 with boundary 
conditions (f(—l) = <f(l) = 0, thus 9 — 1 and 



e k = -g-(fc 

MX) 



sm 



k > 0, 



-(k + l)(X + l) 



(12) 



Since 



cLc x 



(13) 



the particle density of bosons and its correlation function 
© are equal to those of the fermions in the quadratic 
Hamiltonian ([8"]). The fermion two-point function is given 
by 



ab k=0 



4>kx$ky 



exp(w/c/T) 



(14) 



Its TSS limit can be written in terms of the eigensolutions 
of the Schrodinger-like equation (|TU1) , as 



(4 c y) = l ~ 



iv^ (f k (X)(p k (Y) 
^ 1 + exp[(e fc + Hr)/r) ' 



(15) 
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FIG. 2: (Color online) Results for the particle density (top) 
and the density-density correlator (bottom) at fi — 1 and 
p — t- oo, for some values of /, keeping r = Tl 2 = 2 fixed. The 
full line shows the TSS functions given by Eqs. (Tl8|) and ([20) . 



where we have introduced the scaling variable 

t = Tl zB = Tl 26 . (16) 

Then, straightforward calculations allow us to extend the 
T = results of Ref. [25[ , to allow for the temperature. 
We obtain 



p(x) = (n x )=r e V(fil 2e ,Tl 2e ,x/l e ), 
MX) 2 



k=0 



+ exp[(e fc +Ht)/t] 



and 



G(x,y) = r 2e G(fll 2S ,Tl 2S , X /l s ,y/l s ), 

Mx)MY) 



G(jir,T,X,Y) = - 



St 

,fc=0 



+ exp[(efc + Ht)/t] 



(17) 
(18) 

(19) 
(20) 



which holds for x ^ y. Note that G(x, y) < for x ^ y, 
but G{x, x) = (n 2 x ) - (n x ) 2 > 0, indeed Y, xy G(x, y) > 0. 
The above scaling functions describe the asymptotic large 
trap-size behavior; corrections are suppressed by a fur- 
ther I 20 power. Practically exact results for the TSS func- 
tions of the harmonic potential and the hard-wall limit, 
for any fx r and r, can be easily obtained using Eqs. (|TTj) 
and (|12p . because the series are rapidly converging. 



In Fig. Q] we show QMC results for the harmonic po- 
tential at /i = 1, for some values of the trap size, keeping 
t = Tl fixed at r = 2 and r = 8. They clearly appear to 
approach the TSS given by Eqs. (JTHJ) and Analo- 
gous results are obtained for the hard-wall limit, which is 
equivalent to a homogeneous system of size L = 21 with 
open boundary conditions. Since 6 = 1 in this case, TSS 
requires keeping t = Tl 2 fixed. Results for r = 2 are 
shown in Fig. [2j 

The TSS functions in the low-density regime are ex- 
pected to be universal with respect to the on-site coupling 
U, provided that U > 0. As shown by the calculations 
at T = reported in Ref. [30j , finite values of U give rise 
to 0(l~ e ) corrections, which decay more slowly than the 
0{l~ 26 ) leading corrections in the HC limit. 



II. TSS AT THE n = 1 MOTT TRANSITION 

We recall that the quantum critical behavior around 
fi = — 1 of the homogeneous HC BH model without trap 
is essentially analogous to that at p = 1, because of the 
invariance under the particle-hole exchange. At the n = 1 
Mott-insulator to supcrfluid quantum transition, the crit- 
ical exponents v and z and the trap-size exponent 9 are 
the same as those at fj, = 1. However, the particle- hole 
symmetry does not hold in the presence of the trapping 
potential, and the asymptotic trap-size dependence at 
T = appears more complicated at the n = 1 Mott 
transition [25[. 

In the ID HC BH model with a trapping potential 
and \x < 1, the ground-state particle density turns out 
to approach its local density approximation (LDA) in 
the large-? limit, with corrections that are suppressed by 
powers of the trap size and present a nontrivial TSS be- 
haviour [2^], see also below. Within the LDA, the parti- 
cle density at the spatial coordinate x equals the particle 
density of the homogeneous system at the effective chem- 
ical potential 

u eff (x) = fi + {x/lf. (21) 

The LDA of the particle density reads (n x )n a = 
PM*(x/l), where 

!0 for n e s(x) > 1, 

(1/n) arccos/i cff (x) for -1 < p, e s{x) < 1, 
1 for fi cS (x) < -1. 

(22) 

The thermodynamic limit at fixed \x corresponds to the 
limit N, I — >• oo keeping the ratio N/l fixed. At T = we 
have 

N^C£blb t )=p( P )l + 0(l), (23) 

i 

where p(p,) = 2 f °° pida(y) dy. Eq. (|23| provides the cor- 
respondence between the ratio N/l and [i at T — 0. In 
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FIG. 3: (Color online) The rescaled difference of the lowest 
energy levels I 2 A (above) and the rescaled particle density 
I (1— (no)) (below) at the center of the trap vs. / in the middle 
of the trap for p — 2 and at fi — — 1. 

particular, when < N/l < p(— 1) the system is essen- 
tially in the superfluid phase, while for N/l > p(—l) the 
n = 1 Mott phase appears around the center of the trap. 

The T = trap-size dependence around /i = —1 is 
characterized by an infinite number of level crossings of 

(k) 

the two lowest energy states, at Iq with k = 1,2,3,... 

and Zg fc+1 ' > lQ k \ As suggested by Fig. El the rescaled 
energy difference I A of the two lowest states shows a 
periodic structure asymptotically for large /. as a function 
of I [25|]. In particular, in the case of a trap centered at 
the middle site of the lattice and in the large-Z limit, the 
interval 

Pl = l^-l^ k \ (24) 

between two even zeroes of such difference approaches a 
constant value, P* = 1.11072073 for p = 2 and P* = 1 
for p —> 00. This gives rise to a peculiar modulation of the 
amplitudes of the power-law behaviors of the observables. 
For example, the density at the center of the trap behaves 
as 

(n ) - 1 a r%W, (25) 
where 4> is the phase-like variable 

_ l-l Q (2fc ) (2fc+2) , , 

f ~ ,(2fe+2) ,(2fe) ' f - 1 < ( ) l /D J 

'0 'o 

which measures the distance from the closest even level 
crossing. Thus, the amplitudes of the power laws pre- 
dicted by the TSS show a periodic modulation, see 
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FIG. 4: (Color online) The particle density at p = — 1 for 

r = Tl — 8 (top) and r = 2 (bottom). The full lines represent FIG. 5: (Color online) Scaling of the subtracted particle den- 
the LDA approximation (l22l) . sity at p = — 1 for r = 8 (top) and r = 2 (bottom). 



Ref. [25] for details. Moreover, the spatial dependence 
of the particle density and its correlator turns out to be 
described by the following scaling behavior at large trap 
size 



P {x) = (nx) « pMx/i) + r 9 .f(x, 4>), 

where X = x/l 9 , and 

G(0,x) = (n n x ) c = r 2e g(X,<b). 



(27) 



(28) 



An interesting question is whether this scenario persists 
at finite temperature, where the effects of the level cross- 
ings of the lowest states are expected to be weaker. 

We investigate this issue by comparing QMC data at 
various increasing values of the trap size keeping fixed 
the scaling quantity r = Tl. Fig. H] shows results for 
the particle density, which appears to approach the LDA 
approximation (|22[) with increasing I, for r = 2 and r = 
8. This should not be surprising, given that the LDA 
approximation is asymptotically approached at T = 0, 
and the TSS is performed by rescaling T as T — t/1. 
Fig.[5]shows the scaling of the subtracted particle density 



Ap(x) = p{x) - p\da.{x/l). 



(29) 



Its behavior around the trap appears to behave as 

Ap(x) w r e V(X, T ), (30) 
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FIG. 6: (Color online) The density correlation at fj, = — 1 for 
values of the trap size I, at r = Tl = 8 and r — 2. The dashed 
lines show the TSS at p = 1. 



without evidences of modulations. Strong deviations 
from scaling are observed at large distance from the cen- 
ter of the trap, in a region corresponding to x y2l, 
which disappears in the TSS limit, because it moves away 
with increasing / when the data are plotted versus the 
rescaled distance X = x/l 9 with 8 < 1. As we shall see in 
Sec. IIV1 the behavior around x w where p e s{x) w 1, 
shows another scaling behavior, effectively controlled by 
an external linear potential. 
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FIG. 7: (Color online) Plots of J[p(a;) - 1] (top) and l 2 G(0,x) 
(bottom) for p —> oo, at fi — —1 and r = Tl 2 — 2. The 
dashed lines show the TSS at u = 1. 



Some results for the particle-density correlation 
G(0,cc), between the center of the trap and the point 
x, are reported in Fig. [5] They show a nice scaling com- 
patible with 



G(0,x)=r 1 g(X,r), 



(31) 



again without evidence of modulations. Interestingly, the 
data appear to approach the same TSS curves found at 
fi = 1 for the same values of r = TZ, providing an evi- 
dence of universality between the TSS at the n = and 
n = 1 Mott transitions, at least for the connected corre- 
lation functions. Analogous results are found in the limit 
p — > oo, as shown by the data reported in Fig. [7J 

We conclude noting that the apparent absence of mod- 
ulations indicate that the temperature averages out these 
effects, so that the temperature and trap-size depen- 
dences are well described by the TSS ansatz (Q| and ((7} , 
in analogy with the low-density regime. Of course, the 
modulation phenomenon must be somehow recovered in 
the limit r — » 0. 



III. TSS WITHIN THE SUPERFLUID PHASE 

We now analyze the trap-size dependence within the 
japless superfluid phase, whose corresponding continuum 
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FIG. 8: The asymptotic behaviour of the difference between 
the lowest-energy states A and the particle density (no) at 
the center of the trap at fj, = 0, thus / = 1/2, and T — 0. In 
formulae A = -Kt((p)/£ e and (no) = (1 — <fi)/£e, where £ e ~ I" 1 
is the entanglement length scale. These results apply for any 
p, including p — > oo. See Refs. [H, for details. 



theory is a conformal field theory with z = 1. The trap- 
size dependence at zero temperature was investigated in 
Refs. [H, [32| . Its asymptotic behavior turns out to be 
characterized by two length scales with different power- 
law divergence in the large trap-size limit. One of them 
scales as £ ~ / and describes the behavior of observables 
related to smooth modes, such as the half-lattice entan- 
glement; the other one scales as £ ~ V* with ( = p/(p+ 1) 
and it is found in observables involving the modes at the 
Fermi scale kp = nf, where / is the filling of the ho- 
mogeneous system. Moreover, the asymptotic power law 
behaviors appear modulated by periodic functions of the 
trap size, which is again related to level crossings of the 
two lowest states when increasing the trap size, whose 
interval 



Qi = l 



(fc+i) 





,{k) 
'0 i 



(32) 



tends to a constant: Q\ = 1.3110287 for p — 2, and 
Q1 = 1 for p = oo. Indeed, the difference between the 
two lowest states behaves asymptotically as A ~ 
where ip is defined as 



ip = 2 a 



(2k) 



,(2fc+2) _ ,(2fc) ' 



(2k) 



<l <l 



(2k+2) 
i 



(33) 



and t{<p) is a triangle function shown in Fig. [S] where 
results for \i = arc reported. 

At T — 0, the particle density and its correlation func- 
tion show a quite complicated behavior (25|. Indeed, 

(n x ) w pu,{X) + r^hiY, V )e 2ikFX + g(Y, <p)}, 
X = x/l, Y = xl~ p/{p+1 \ (34) 



and 



G 



n (x,0) « r 2 *l<*+V%t\h{Y,<p)e 2ik * x +g(Y,<p)}, (35) 
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FIG. 9: (Color online) We plot the particle density (top), the 
rescaled subtracted particle density lAp with Ap = p — pi da 
(middle), and the particle density correlator G(0, x) (bottom) 
at p — and T = 2/1, in the case of the harmonic potential. 
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FIG. 10: (Color online) Scaling of the particle density (top), 
the subtracted particle density Ap = p—pida, (middle) and the 
particle density correlator (bottom) at p — and Tl 2 ^ 3 = 1 
around x — I where p c g = 1. 



where kp = nf = arccos/j, and terms suppressed by 
higher powers of Z _1 are neglected. 

In Fig. [9] we show results for the particle density and 
various values of the trap size with T = 2/1. We choose 
this scaling of the temperature because the trap expo- 
nent associated with the smooth modes is expected to 
be 9 = 1 [23. The results are clearly converging toward 
the corresponding LDA approximation. Around the cen- 
ter of the trap, we do not observe the modulated scaling 
found at T = 0. Only around x ~ I we observe some sig- 
nificant differences which decrease with increasing I, see 
below. Moreover, the data for the particle-density corre- 
lator G(0, x) appear to vanish after a few lattice spacings 
for any I, without showing any particular scaling. There- 



fore, the leading T = scaling behavior (|35|) associated 
with the exponent 2p/(p+ 1) is also effectively averaged 
out by the temperature. 

IV. TSS AROUND THE SPATIAL POINT 
WHERE p cH {x) = 1 

When fi < 1 , the region where the particle density ap- 
pears to rapidly vanish corresponds to the region where 
the effective chemical potential defined in Eq. (|21l) ap- 
proaches the values p, e g ~ 1, which is the value of the 
chemical potential corresponding to the transition be- 
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tween the superfluid phase and the empty state, where 
the particle density of the ground state gets suppressed. 
We thus expect that, for generic values of p and p, the 
region around 



x c = 1(1 - p) 



i/p 



(36) 



where p e g(x c ) — 1, develops critical modes related to the 
superfluid to vacuum transition. The effective chemical 
potential can be expanded around 

/irfOr) = 1+ P (l- ^(P-D/f £^ + 0[(x - x c f]. (37) 

Therefore the behavior around x c is essentially analogous 
to that arising at p = 1 in the presence of a linear poten- 
tial Vi ~ r/l. Around x c , critical modes should appear 
with length scale £ ~ Z CT , where a is the exponent associ- 
ated with a linear external potential, thus [25[ a = 1/3, 
obtained by replacing p = 1 in Eq. ([5} . We then expect 
that around x — x c 

P (x) = l-Wfplix - z c )/i 1/s , TZ 2 / 3 ], (38) 
Ap(x) =p- p lda = r 1/3 f Ap [(x - x,)// 1 / 3 , Tf/% (39) 

(40) 



G(x c , x) = l-Wfg[(x - as c )/i 1/8 , TZ 2/3 



In Fig.[TU]we show QMC results at p = for some values 
of Z keeping TZ 2 / 3 = 1 fixed, which nicely confirm the 
scaling behaviors (138H40[) . Note that analogous results 
are obtained for any p < 1, around the region where 

p e s(x) « 1. 



V. CONCLUSIONS 

We have investigated the interplay of temperature and 
trap effects in cold particle systems at their quantum crit- 
ical regime, such as cold bosonic atoms in optical lattices 
at the transitions between Mott-insulator and superfluid 
phases. The theoretical framework is provided by the ID 
Bose-Hubbard model ([1} in the presence of an external 
trapping potential, and by the TSS theory leading to the 
TSS ansatz Q and J7]). This study is of experimental rel- 
evance, because quasi ID trapped particle systems have 
been realized in optical lattices, see, e.g., Refs. [Jill. 
Therefore, temperature and trap effects at the Mott tran- 
sitions may be investigated by improving the accuracy of 
the experiments. 

Our numerical study is based on QMC simulations 
at fixed chemical potential and trap size. Since the 
main features of the TSS at the Mott transitions are ex- 
pected to be universal with respect to the on-site cou- 
pling U (provided U > 0), we consider the HC limit 
U — > oo, which is particularly convenient because it re- 
stricts the values of the site particle number to n x = 0, 1, 
and it is expected to minimize the corrections to the 
universal TSS, as verified at T — [3(| • We present 
finite-temperature results for the particle density and the 



density-density correlator, at the quantum Mott transi- 
tions, and within the gapless superfluid phase. Their 
temperature and trap-size dependences appear well de- 
scribed by the simplest TSS ansatz, such as that given 
in Eq. ([7]) with 8 = p/(p + 2). In particular, the finite- 
temperature data do not show evidence of the periodic 
modulations of the asymptotic behaviors which are ob- 
served at T = [25| . They are effectively averaged out 
by a nonzero temperature, leaving only a plain power-law 
TSS behavior. 



Acknowledgments 

The QMC simulations were performed at the INFN 
Pisa GRID DATA center, using also the cluster CSN4. 



Appendix A: The QMC simulations 

Our numerical study of the BH model (JTJ) is based on 
QMC simulations using the stochastic series expansion 



method [27[ with the directed loop algorithm [28[. We 
compute the particle density (n x ) and its connected cor- 
relation function (n x n y ) c . 

The QMC simulations are performed at fixed tempera- 
ture T, chemical potential p, trap size Z, and lattice size L 
with open boundary conditions. In order to avoid finite- 
size effects, we choose L sufficiently large to effectively 
obtain L^oo data within the statistical errors. 

QMC simulations require also to fix some algorithm 
parameters, such as the order A tr of the truncation of 
the Taylor expansion of the partition function, and the 
number of operator loops N\ oops . We set N tr in the stan- 
dard way as A tr = 1.5 M max N honds /T where M max is the 
highest matrix element of the single bond Hamiltonians 
and Abonds is the number of interacting site pairs. At 
the end of each simulation we analyze the fluctuation in 
the order of the series expansion, and verify whether the 
averaged order is safely less than A tr , and its deviations 
have the expected behavior proportional to the square 
root of the mean value. We define a MC step as a single 
sweep of identity and diagonal operators in the expansion 
of the partition function followed by a number Ai oops of 
directed operator loop updates. The number Ai oops is 
fixed by requiring that the number of visited vertices is 
globally of the order of A tr in a MC step. 

Typical statistics of our QMC simulations are a few 
million MC steps for each set of parameters. We have 
checked that, in all QMC simulations presented here, the 
autocorrelation times of the observables that we compute 
are smaller than 10 2 in units of MC steps. Effective L — > 
oo data within the resulting statistical errors are obtained 
by taking L/l as 4 for the simulations at p — 1 and p = 0, 
and L/l ss 9 for p = — 1. 

Let us now discuss the operator- loop equations. Since 
the Hamiltonian of the HC BH model can be written 
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as a sum of bond terms, and since there is a one-to- 
one correspondence between the site occupation number 
and the states of the XXZ model [HI, the non-zero ma- 
trix elements are the same as those reported in Ref. [28[ 
(see its Sec. 2). Of course, their values are different, in 
particular the diagonal matrix elements are site depen- 
dent, due to the existence of the trap which gives rise 
to a space-dependent effective chemical potential fi c g(x). 
Using considerations based on the symmetry and con- 
servation laws, we can divide the transition probabilities 



involved in the operator-loop update in eight indepen- 
dent groups for each site, as in Ref. [28| (see its Sec. 3), 
with three possible bounces (probabilities to exit the ac- 
tual visited vertex of the operator loop from the same 
entrance leg) for each group. The corresponding system 
of equations are solved by minimizing the bounces. As 
a result, we allow one bounce within each group in the 
QMC simulations at fj, = 1 and \x = 0, and two bounces 
at fi = —1. 
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